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ABSTRACT: Radio astronomy has provided evidence for the presence of ionized atmospheres around almost all 
classes of non-degenerate stars. Magnetically confined coronae dominate in the cool half of the Hertzsprung-Russell 
diagram. Their radio emission is predominantly of non-thermal origin and has been identified as gyrosynchrotron 
radiation from mildly relativistic electrons, apart from some coherent emission mechanisms. Ionized winds are 
found in hot stars and in red giants. They are detected through their thermal, optically thick radiation, but 
synchrotron emission has been found in many systems as well. The latter is emitted presumably by shock- 
accelerated electrons in weak magnetic fields in the outer wind regions. Radio emission is also frequently detected 
in pre-main sequence stars and protostars, and has recently been discovered in brown dwarfs. This review 
summarizes the radio view of the atmospheres of non-degenerate stars, focusing on energy release physics in cool 
coronal stars, wind phenomenology in hot stars and cool giants, and emission observed from young and forming 
stars. 

Eines habe ich in einem langen Leben gelernt, namlich, dass unsere 
ganze Wissenschaft, an den Dingen gemessen, von kindlicher Primi- 
tivitdt ist - und doch ist es das Kostlichste, was wir haben. 
One thing I have learned in a long life: that all our science, measured 
against reality, is primitive and childlike - and yet it is the most precious 
thing we have. 

A. Einstein 1951, in a letter to H. Miihsam, Einstein Archive 36-610 



1 INTRODUCTION 



Stellar radio astronomy has matured over the past two decades, driven in particular by discoveries 
made with the largest and most sensitive radio interferometers. Radio emission is of great 
diagnostic value as it contains telltale signatures not available from any other wavelength regime. 
Some of the detected radio emission represents the highest-energy particle populations (MeV 
electrons) yet accessible on stars, the shortest (sub-second) detectable time scales of variability 
and energy release, and probably refers most closely to the primary energy release responsible 
for coronal heating. This review is to a large extent devoted to demonstrating the ubiquity of 
high-energy processes in stars as revealed by radio diagnostics. 

Stellar radio sources include thermal and non-thermal magnetic coronae, transition regions 
and chromospheres, stars shedding winds, colliding-wind binaries, pre-main sequence stars with 
disks and radio jets, and embedded young objects visible almost exclusively by their radio and 
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millimeter-wave emission. Most recent additions to the zoo of objects are brown dwarfs, and 
with the increasingly blurred transition from stars to brown dwarfs to giant planets like Jupiter 
and Saturn, even the magnetospheres of the latter may have to be considered a manifestation 
of magnetic activity in the widest sense. To keep the discussion somewhat focused, this paper 
concentrates on physical processes in magnetic coronae, but includes, in a more cursory way, 
atmospheres of young and forming stars and winds of hot stars. We do not address in detail 
the large phenomenology of extended and outflow-related sources such as radio jets, HII regions, 
masers, Herbig-Haro objects, and the diverse millimeter /submillimeter phenomenology, e.g., 
molecular outflows and dust disks. Compact stellar objects (white dwarfs, neutron stars, black 
holes) are not considered here. 

Inevitably - and fortunately - much of the knowledge gained in stellar astronomy is anchored 
in solar experience. The privilege of having a fine specimen - and even an exemplary prototype 
- next door is unique among various fields of extrasolar astrophysics, being shared since recently 
only by the related field of extrasolar planetary astronomy. Detailed solar studies, even in-situ 
measurements of the solar wind, are setting high standards for studies of stellar atmospheres, 
with a high potential reward. Solar radio astronomy has been reviewed extensively in the 
literature. For detailed presentations, we refer to Dulk (1985) and Bastian et al. (1998). 

The maturity of stellar radio astronomy is demonstrated by a number of review articles on 
various subjects; a non-exhaustive list for further reference includes Andre (1996), Bastian (1990, 
1996), Bookbinder (1988, 1991), Dulk (1985, 1987), Giidel (1994), Hjellming (1988), Kuijpers 
(1985, 1989a), Lang (1990, 1994), Lestrade (1997), Linsky (1996), Melrose (1987), Mullan (1985, 
1989), Mutel (1996), Phillips (1991), Seaquist (1993), van den Oord (1996), and White (1996, 
2000). A natural starting point for this review, even if not strictly adhered to, is Dulk's com- 
prehensive 1985 Annual Review article that summarizes the pre- and early- Very Large Array 
view of stellar (and solar) radio emission. Meanwhile, two dedicated conferences, the first one 
in Boulder in 1984 (proceedings edited by Hjellming & Gibson 1985) and the second held in 
Barcelona in 1995 (proceedings edited by Taylor & Paredes 1996), provided a rich forum to 
discuss new developments; together, they beautifully illustrate the progress made over the past 
decades. 

2 RADIO SURVEYS AND THE RADIO HERTZSPRUNG-RUSSELL DI- 
AGRAM 

Figure 1 presents a radio Hertzsprung-Russell Diagram (HRD) based on stellar radio detections 
between 1—10 GHz, reported in the catalog of Wendker (1995); for other examples, see White 
(1996, 2000). The luminosity is the logarithmic average of all reported detections. The accuracy 
of the location of some stars on the HRD is compromised by the limited quality of distance and 
color measurements, or in multiple systems, by the uncertain attribution of the radio emission to 
one of the components. Nevertheless, almost all of the usual features of an HRD are recovered, 
testifying to the importance and ubiquity of radio emission. 

The many nearby M dwarfs in Figure 1 were among the first radio stars surveyed. The 
typically larger distances to earlier-type cool dwarfs made their discovery more of an adventure, 
but the samples now include late-type binaries, K stars, solar analog G stars, and a few F stars. 
Most radio-detected dwarf stars are strong X-ray sources and young, rapid rotators. 

The cool half of the subgiant and giant area is dominated by the large and radio-strong sample 
of RS CVn and Algol-type close binaries. Other sources in this area include the vigorous coronal 
radio sources of the FK Com class, chromospheric radio sources, symbiotic stars and thermal- 
wind emitters. Additionally, a very prominent population of thermal or non-thermal sources 
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Figure 1: HR diagram showing 440 radio-detected stars. 



just above the main sequence is made up of various classes of pre-main sequence stars, such as 
classical and weak-lined T Tauri stars. 

Common to most cool star radio emitters is their non-thermal nature which, especially in its 
"quiescent" form, constitutes one of the most significant - and unexpected - discoveries in stellar 
radio astronomy, suggesting the presence of magnetic fields and high-energy electrons. Moving 
toward A-type stars on the HRD, one expects, and finds, a dearth of radio detections owing 
to the absence of magnetic dynamo action. However, this is also the region of the chemically 
peculiar Ap /Bp stars that possess strong magnetic fields and many of which are now known to 
be non-thermal radio sources as well. Some pre-main sequence Herbig Ae/Be stars in this area 
are also prominent radio emitters. 

Finally, the hot-star region is heavily populated by luminous wind-shedding O and B stars 
and Wolf-Rayet (WR) stars, both classes showing evidence for either thermal or non-thermal 
radio emission. 
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3 THEORY OF RADIO EMISSION FROM STELLAR ATMOSPHERES 

3. 1 Elementary Formulae for Radiation and Particles 

We summarize below some handy formulae that have become standard for most radio-stellar 
interpretational work (the interested reader is advised to consult the original literature as 
well). In thermodynamic equilibrium, the emissivity of a plasma of temperature T e gf, rj v (in 
erg s~ 1 cm _3 Hz _1 sr~ 1 ), and the absorption coefficient, k u (in cm' 1 ) are related by Kirchhoff's 
law; for combined modes of polarization, for unity spectral index and a frequency v <C 10 10 T e g 
(y in Hz, T c g in K; Rayleigh- Jeans approximation): 

* = ~ 3.1 x 10^T eS u 2 . (1) 

K v C 

Here, k = 1.38 x 1CT 16 erg K _1 is the Boltzmann constant. The brightness temperature T& is 

rp J t T cff , r < 1 , , 

Tb = \T cS , r»l ® 

with t = J ndl being the optical depth along the line of sight t. The spectral radio flux density 
S u then is 

^ = ^^Ki7^)(dk) (iot^) (cf) mJy ^ 

where A is the cross sectional source area perpendicular to the line of sight, and the approxima- 
tion on the right-hand side is for a circular source with A = Trr 2 . In non-thermal astrophysical 
plasmas, electron volumetric number densities are often observed to follow a power-law form 

n(e) = N(S - l^V 5 [cm^crg- 1 ] (4) 

where e = (7 — l)m e c 2 is the kinetic particle energy, 7 is the Lorentz factor, and 5 > 1 has been 
assumed so that iV is the total non-thermal electron number density above eo- A fundamental 
frequency of a plasma is its plasma frequency, given by 

o2\ V2 



Up = 



^P_(n^ ] ^gooony 2 [Hz]. (5) 



2tt \ 7rm f 



Here, m e = 9.1 x 10 28 g is the electron rest mass, e = 4.8 x 10 10 esu is the electron charge, 
and n e is the total electron number density (cm -3 ). 



3. 2 Bremsstrahlung 

A thermal plasma emits free-free emission (bremsstrahlung) across the electromagnetic spectrum. 
For cosmic abundances, the absorption coefficient is then approximately given by (in the usual 
units; T = T c g, mix of 90% H and 10% He; upper equation for singly ionized species, lower 
equation for fully ionized plasma, and v 3> v v \ Dulk 1985, Lang 1999): 

k ~ 01n 2 T- 3 / 2 ^ 2 x / ln ( 5 X 107T3/ >), ^ ^ 3.2 x 10 5 K ( ) 

- u - Ui ^M v x I ln(4 7 x 10 io r/l/)> T ^ 3.2 x 10 5 K [b) 

and the emissivity n u follows from equation 1. Optically thick bremsstrahlung shows a u 2 
dependence, whereas optically thin flux is nearly independent of v. A homogeneous optically 
thin magnetized source is polarized in the sense of the magnetoionic x-mode, whereas an optically 
thick source is unpolarized (Dulk 1985). 
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3.3 Gyromagnetic Emission 

Electrons in magnetic fields radiate gyromagnetic emission. The gyrofrequency in a magnetic 
field of strength B is 

u c = — = « 2.8 x 10 6 B [Hz] (7) 

where c = 3 x 10 10 cm s" 1 is the speed of light, and the magnetic field strength B is given in 
Gauss. The relativistic gyrofrequency is f Circ i = v c /l- For large pitch angles, spectral power is 
predominantly emitted around a harmonic s with 

fmax = S^c.rcl ~ T^rel ~ 2.8 X 10 6 J B 7 2 [Hz] . (8) 

Depending on 7, the emission is termed gyroresonance or cyclotron emission (s < 10, 7 « 1, 
non-relativistic, typically thermal electrons), gyro synchrotron (s « 10 — 100, 7^2 — 3, mildly 
relativistic electrons), or synchrotron emission (s > 100, 7> 1, relativistic electrons). Because 
relativistic effects play an increasingly important role in increasing s, the fundamental properties 
(e.g., directivity, bandwidth, polarization) change between the three categories of gyromagnetic 
emission. The total power emitted by an electron is (Lang 1999) 

P = 1.6 x 10" 15 /3 2 7 2 B 2 sin 2 a [erg s" 1 ] (9) 

where (3 = v/c = (1 — 7 -2 ) 1 / 2 and a is the pitch angle of the electron. 

Simplified approximate expressions for separate magnetoionic modes of gyromagnetic emission 
have been given by Dulk & Marsh (1982), Robinson & Melrose (1984), and Klein (1987). For 
stellar observations, the angle 6 between the line of sight and the magnetic field is often unknown 
and should be averaged. We give handy, simplified expressions for n u , k v , the turnover frequency 
^peak, and the degree of circular polarization r c derived from Dulk (1985) for = 7r/3, with 
exponents for the model parameters B, N, T, and the characteristic source scale along the line 
of sight, L. For more comprehensive expressions, we refer to the original work by Dulk & Marsh 
(1982) and Dulk (1985). 

Gyroresonance emission: If s is the harmonic number, s 2 T «6x 10 9 K and v ^ u c , then for 
each magnetoionic mode (for unity refractive index), 



k v (s) = 1020(1 ±0.5) 2 n& * 2 I "~ l 



vT 1 / 2 s\ \ 1.6 x 10 10 




(1 - svjuf 



8.4 x 10 



-li' 



(10) 



Vu(s) = ^-^-Ku(s) (11) 

(for the x-mode, and marginally for the o-mode, see Dulk 1985) where n e is the ambient elec- 
tron density. Properties: Emission is strongly concentrated in emission "lines" at harmonic 
frequencies sf c , and k u depends on the magnetoionic mode (lower sign for o-mode, upper sign 
for x-mode). The emission at a given harmonic comes from a layer of constant B with a 
thickness L = 2Ab/3o c os# ~ Ab/?o determined by the magnetic scale length Ab = B/VB and 
fa = [kT/{m e c 2 )} 1 / 2 w 1.3 x 10- 5 TV2. 
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Gyro synchrotron Emission from a Thermal Plasma: If 10 s ^ T ^ 10 9 K and 10 ^ s ^ 100, then 
for the x-mode (and r c for r <C 1) 



K,, 



i 21n e T 7 B 9 u- 10 (12) 

« 3.2 x 10" 36 n e T 8 SV" 8 (13) 

*Vak ~ 1.3{n e L/B) 01 T°- 7 B (14) 

: 2.9 x 10 4 T- ai3 V/5)- a51 . (15) 



r. 



Properties: The spectral power is oc v 2 on the optically thick side, but oc v 8 on the optically 
thin side of the spectrum for a homogeneous source. 

Gyro synchrotron Emission from a Power-Law Electron Distribution: For isotropic pitch angle 
electron distributions according to equation 4, with eo = 10 keV = 1.6 x 10~ 8 erg, 2 <^ 5 & 7, 
and 10 & s & 100, and for the x-mode (and r c for r«l) 



K„ « 10 -0.47+6.0M iVjB 0.3+0.98* I/ -1.3-0.9M (1 g) 

^ » 10 -31.32+5.2« JVB -0.22+0.9« l/ 1.22-0.9tf (17) 

^pcak ~ 1 3-41+0.275 (7VL) 0.32-0.035 B 0.68+0.035 (lg) 

r c « 10 3.35 + 0.0355 (z//s) -0.51 

Properties: Broad spectra with intermediate circular polarization degree on the optically thin 
side. Optically thick spectral power is approximately oc z/ 5 / 2 , and the optically thin side is a 
power law with an index a = 1.22 — 0.95 for a homogeneous source. 

Synchrotron Emission from a Power-Law Electron Distribution (homogeneous and isotropic): 
For 7> 1, i.e., s S> 100, in each of the magnetoionic modes 

K„ « 10 5.89+1.7M ((J _ 1)iVB (*+2)/2 I/ -(*+4)/2 (2Q) 
Tfc, « 10 -24.7+1.575 ((5 _ 1)ivB (5 + l)/2^-(5-l)/2 (21) 

^peak '- 



10 11.77+3.44<5 (J _ 1 )2^2 L 2 B *+2l V(*-M) (22) 



Properties: Continuous and broad spectrum, important harmonics s « 7 3 . For a homogeneous, 
optically thin source, the degree of linear polarization is r = (<5 + l)/(5 + 7/3). The optically 
thick (self-absorbed) spectral power is oc z^ 5 / 2 , and the optically thin side of the spectrum has a 
power-law spectral index of a = —(5 — l)/2. 



3.4 Coherent Emission 

The brightness temperature of synchrotron emission is limited to T b £ 10 12 K by inverse Comp- 
ton scattering (Kellermann & Pauliny-Toth 1969). If higher T5 is observed, then a coherent 
radiation mechanism should be considered. Two mechanisms have received most attention both 
for solar and stellar coherent bursts (for details, see Benz 2002): 

Plasma radiation is emitted at the fundamental or the second harmonic of v v (equation 5; for 
propagation parallel to the magnetic field), or of the upper hybrid frequency is\jh = {yl, + v 2 ) 1 / 2 
(for perpendicular propagation, with v p > u c ). It is thus a useful means to approximately de- 
termine the electron density in the source. It can account for high brightness temperatures 
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(up to 10 18 K, Melrose 1989) and small bandwidth, and it is frequently observed in the Sun 
(Bastian et al. 1998) and in low- frequency stellar flares (Bastian 1990). Owing to increasing 
free-free absorption with increasing v (because n e oc in equation 6), and possibly also owing 
to gyroresonance absorption in hot plasma, fundamental emission is best observed below 1 GHz, 
although in very-high temperature environments such as coronae of RS CVn binaries, the ab- 
sorption is milder and the limitations are more relaxed (White & Franciosini 1995). The escape 
of fundamental emission is also alleviated in a highly structured medium with sharp gradients 
in density (Aschwanden 1987, Benz et al. 1992). 

Electron cyclotron maser emission (Melrose & Dulk 1982, 1984) is emitted mostly at the 
fundamental and the second harmonic of v c (equation 7) and can therefore be used to determine 
the magnetic field strength in the source. The requirement for radiation propagation, sv c > u p , 
also sets an upper limit to the electron density in the source and along the line of sight to the 
observer. The cyclotron maser mechanism accounts for the observed high TJ, ( ^ 10 20 K for 
s = 1, ^ 10 16 ~ 17 K for s = 2), and polarization degrees up to 100%. 

Other mechanisms have occasionally been proposed. For example, plasma maser emission 
from a collection of double layers was initially calculated by Kuijpers (1989b) and applied to 
highly polarized stellar flare emission by van den Oord & de Bruyn (1994). 



3.5 Wind emission 

Olnon (1975), Panagia & Felli (1975), and Wright & Barlow (1975) calculated the radio spectrum 
of a spherically symmetric, ionized wind. If M is the mass loss rate in units of M yr _1 , v the 
terminal wind velocity in km s _1 , d pc the distance to the star in pc, and T the wind temperature 
in K, the optically thick radius is (we assume a mean atomic weight per electron of 1.2, and an 
average ionic charge of 1) 

^thick = 8xl0 28 (M/t;) 2 / 3 T- a45 ^ - 7 [cm]. (23) 
Then, the following formulae apply (R* in cm): 

/ 9 x lO^M/^/ST - 1 ^ - 6 My], if ^thick > R* (0A , 

V [ 5 x 10 39 (MA) 2 T-°- 35 i^p>- ai [mJy], if iW < R*. 1 ] 



3.6 Loss Times 

Relativistic electrons in an ambient gas of density n e lose energy by Coulomb collisions with the 
ions. The energy loss rate and the corresponding life time are (under typical coronal conditions; 
Petrosian 1985) 

-Tcoii = 5 x 10" 13 n e [s" 1 ], 7-cou = 2 x 10 12 ^- [s]. (25) 

n e 

For electrons in magnetic fields, the synchrotron loss rate and the life time are given by (Petrosian 
1985, pitch angle = vr/3) 

-7B = 1.5 x 10- 9 J5V Is" 1 ], t b = ^-^- [s]. (26) 
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4 RADIO FLARES FROM COOL STARS 

4-1 Incoherent Flares 

Many radio stars have been found to be flaring sources, and a considerable bibliography on 
stellar radio flares is now available. As in the Sun, two principal flavors are present: incoherent 
and coherent radio flares. Incoherent flares with time scales of minutes to hours, broad-band 
spectra, and moderate degrees of polarization are thought to be the stellar equivalents of solar 
microwave bursts. Like the latter, they show evidence for the presence of mildly relativistic 
electrons radiating gyrosynchrotron emission in magnetic fields. Many flares on single F/G/K 
stars are of this type (Vilhu et al. 1993, Giidel et al. 1995b, 1998), as are almost all radio flares 
on M dwarfs above 5 GHz (Bastian 1990, Giidel et al. 1996a), on RS CVn binaries (Feldman 
et al. 1978, Mutel et al. 1985), on contact binaries (Vilhu et al. 1988), and on other active 
subgiants and giants although some of the latter perhaps stretch the "solar analogy" beyond 
the acceptable limit: The FK Corn-type giant HD 32918 produced flare episodes lasting 2—3 
weeks, with a radio luminosity of 6 x 10 19 erg s _1 Hz _1 (Slee et al. 1987b, Bunton et al. 1989). 
Its integrated microwave luminosity is thus about 1000 times larger than the total X-ray output 
of the non-flaring Sun! 

4-2 Coherent Radio Bursts 

Bursts that exhibit the typical characteristics of coherent emission (see "Coherent Emission," 
above) probably represent stellar equivalents of metric and decimetric solar bursts that them- 
selves come in a complex variety (Bastian et al. 1998). Some exceptionally long-lasting (~1 h) 
but highly polarized bursts require a coherent mechanism as well (examples in Lang & Willson 
1986a, White et al. 1986, Kundu et al. 1987, Willson et al. 1988, and van den Oord & de Bruyn 
1994). Coherent flares are abundant on M dwarfs at longer wavelengths (20 cm; Kundu et al. 
1988, Jackson et al. 1989). There are also some interesting reports on highly polarized, possibly 
coherent bursts in RS CVn binaries (Mutel & Weisberg 1978, Brown & Crane 1978, Simon et al. 
1985, Lestrade et al. 1988, White & Franciosini 1995, Jones et al. 1996, Dempsey et al. 2001). 
After many early reports on giant metric flares observed with single-dish telescopes (Bastian 
1990), the metric wavelength range has subsequently been largely neglected, with a few notable 
exceptions (Kundu & Shevgaonkar 1988, van den Oord & de Bruyn 1994). 

Coherent bursts carry profound information in high-time resolution light curves. Radio "spike" 
rise times as short as 5—20 ms have been reported, implying source sizes of r < cAt s=s 1500- 
6000 km and brightness temperatures up to ~ 10 16 K, a clear proof of the presence of a 
coherent mechanism (Lang et al. 1983, Lang & Willson 1986b, Giidel et al. 1989a, Bastian et 
al. 1990). Quasi-periodic oscillations were found with time scales of 32 ms to 56 s (Gary et al. 
1982, Lang & Willson 1986b, Bastian et al. 1990), and up to «5 min during a very strong flare 
(Brown & Crane 1978, although this emission was proposed to be gyrosynchrotron radiation). 

4-3 Radio Dynamic Spectra 

If the elementary frequencies relevant for coherent processes, v v and u c , evolve in the source, or 
the radiating source itself travels across density or magnetic field gradients, the emission leaves 
characteristic traces on the u — t plane, i.e., on dynamic spectra (Bastian et al. 1998). A rich 
phenomenology has been uncovered, including: a) short, highly polarized bursts with structures 
as narrow as Avjv = 0.2% suggesting either plasma emission from a source of size ~ 3 x 10 8 cm, 
or a cyclotron maser in magnetic fields of ~ 250 G (Bastian k. Bookbinder 1987, Giidel et al. 
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Figure 2: Gallery of radio dynamic spectra of M dwarf flares. Upper three panels show a flare 
on AD Leo, recorded with the Arecibo (top), Effelsberg (middle) and Jodrell Bank (bottom) 
telescopes in different wavelength ranges (see also Giidel et al. 1989a). Bottom three panels show 
flares on AD Leo (top and middle) and YZ CMi (bottom) observed at Arecibo (after Bastian et 
al. 1990, reproduced with permission of the AAS). 



1989a, Bastian et al. 1990). b) Evidence for spectral structure with positive drift rates of 2 
MHz s _1 around 20 cm wavelength, taken as evidence for a disturbance propagating "down- 
ward" in the plasma emission interpretation (Jackson et al. 1987); and c) in solar terminology, 
rapid broadband pulsations, "sudden (flux) reductions", and positive and negative drift rates of 
250-1000 MHz s" 1 (Bastian et al. 1990, Abada-Simon et al. 1994, 1997, see Figure 2). 

The smallest spectral bandwidths were found to be in the 1% range for some bursts. Conser- 
vatively assuming a magnetic scale height of A# = li?*, the source size can be estimated to be 
r ~ {/S.v /v)A.b ~ a few 1000 km, again implying very high Tj, (Lang & Willson 1988a, Giidel et 
al. 1989a, Bastian et al. 1990). 
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5 QUIESCENT EMISSION FROM CORONAE OF COOL STARS 

5. 1 Phenomenology 

The discovery of quiescent radio emission from magnetically active stars was one of the most 
fundamental and one of the least anticipated achievements in radio astronomy of cool stars 
because there simply is no solar counterpart. The Sun emits steady, full-disk, optically thick 
thermal radio emission at chromospheric and transition region levels of a few 10 4 K. However, one 
derives from equation 3 that such emission cannot be detected with present-day facilities, except 
for radiation from the very nearest stars, or giants subtending a large solid angle (see "Radio 
Emission from Chromospheres and Winds," below). The radio luminosity of the corona caused 
by optically thin bremsstrahlung is proportional to the observed X-ray emission measure, but 
again the calculated radio fluxes are orders of magnitude below observed fluxes (Gary & Linsky 
1981, Topka & Marsh 1982, Kuijpers & van der Hulst 1985, Borghi k Chiuderi Drago 1985, Massi 
et al. 1988, Giidel & Benz 1989, van den Oord et al. 1989). Obviously, another mechanism is in 
charge, and its characterization requires a proper description of the phenomenology. 

Quiescent emission is best defined by the absence of impulsive, rapidly variable flare-like 
events. Common characteristics of quiescent emission are (i) slow variations on time scales 
of hours and days (Pallavicini et al. 1985, Lang & Willson 1986a, 1988a, Willson et al. 1988, 
Jackson et al. 1989), (ii) broad-band microwave spectra (Mutel et al. 1987, Giidel & Benz 1989), 
(hi) brightness temperatures in excess of coronal temperatures measured in X-rays (Mutel et al. 
1985, Lang k Willson 1988a, White et al. 1989a), and often (iv) low polarization degrees (Kundu 
et al. 1988, Jackson et al. 1989, Drake et al. 1992). Occasionally, strong steady polarization up 
to 50% (Linsky k Gary 1983, Pallavicini et al. 1985, Willson et al. 1988, Jackson et al. 1989) or 
unexpectedly narrow-band steady emission (Lang k Willson 1986a, 1988a) is seen on M dwarfs. 
Quiescent emission has been reported at frequencies as low as 843 MHz (Vaughan k Large 1986, 
Large et al. 1989) and as high as 40 GHz (White et al. 1990). 

5.2 Gyromagnetic Emission Mechanisms 

Because active stars show high coronal temperatures and large magnetic filling factors (Saar 
1990) that prevent magnetic fields from strongly diverging with increasing height, the radio 
optical depth can become significant at coronal levels owing to gyroresonance absorption (equa- 
tion 10). Such emission is also observed above solar sunspots. The emission occurs at low 
harmonics of u c , typically at s = 3 — 5 (Giidel k Benz 1989, van den Oord et al. 1989, White et 
al. 1994). Rising spectra between 5 and 22 GHz occasionally observed on stars support such an 
interpretation (Cox k Gibson 1985) and imply the presence of ~ 10 7 K plasma in strong, kG 
magnetic fields (Giidel k Benz 1989; Figure 3). However, in most cases high-frequency fluxes are 
lower than predicted from a uniform, full-disk gyroresonance layer, constraining the filling factor 
of strong magnetic fields immersed in hot plasma. It is possible that the hot plasma resides in 
low-B field regions where the gyroresonance absorption is negligible, for example between strong 
fields from underlying active regions, or at large heights (White et al. 1994, Lim et al. 1996a). 
Uniform magnetic structures of ~600 G containing hot plasma can, on average, not reach out 
to beyond 1-2.R* on M dwarfs (Leto et al. 2000). 

The gyroresonance mechanism cannot apply to lower-frequency radio emission, because the 
radius of the optically thick layer, still at s = 3 — 5, would have to be > 3R* for dMe stars. 
However, an extrapolation of the corresponding magnetic fields of more than 100 G down to 
photospheric levels would result in photospheric fields much stronger than observed (Gary k 
Linsky 1981, Topka k Marsh 1982, Linsky k Gary 1983, Pallavicini et al. 1985, Lang k Willson 
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Figure 3: Radio spectra of the RS CVn binary HR 1099 (upper set) and of the dMe dwarf UV 
Cet (lower set) at different flux levels. The gently bent spectra are indicative of gyrosynchrotron 
emission, and the high-frequency part of U-shaped spectra for UV Cet has been interpreted as 
a gyroresonance component (HR 1099 spectra: courtesy of SM White). 

1986a, Willson et al. 1988, Giidel & Benz 1989). 

The only remedy is to allow for much higher T e g. Then, the optically thick layer shifts to 
harmonics above 10, the range of gyrosynchrotron radiation, and the optically thick source 
sizes become more reasonable for M dwarfs, R ~ i?* (Linsky & Gary 1983, Drake et al. 1989, 
1992). For a thermal plasma, however, the spectral power drops like u~ s at high frequencies 
(equation 13), far from the rather shallow ^-(O- 3 --- 1 ) spectra of active stars (Massi & Chiuderi- 
Drago 1992; Figure 3). On the other hand, an optically thick thermal contribution requires strong 
magnetic fields (~200 G) over such large source areas that inferred photospheric magnetic fields 
again become unrealistically large (Kuijpers & van der Hulst 1985, Drake et al. 1992). 

The situation is much more favorable for a non-thermal electron energy distribution such as 
a power law, analogous to distributions inferred for solar and stellar microwave flares (Kundu & 
Shevgaonkar 1985, Pallavicini et al. 1985). This model is supported by measured high brightness 
temperatures (Lestrade et al. 1984b, Mutel et al. 1985, Umana et al. 1991, Benz et al. 1995). 
Comprehensive spectral modeling suggests mildly relativistic electrons in ~100 G fields with 
power-law indices of 5 « 2 — 4, matching observed broad spectra with turnover frequencies in 
the 1 - 10 GHz range (White et al. 1989b, Slee et al. 1987a, 1988, Umana et al. 1998; Figure 3). 
The non-thermal model is now quite well established for many classes of active stars and provokes 
the question of how these coronae are continuously replenished with high-energy electrons. 

Occasionally, steady narrow-band or strongly polarized coherent emission has been observed at 
low frequencies. An interesting possibility is large numbers of unresolved, superimposed coherent 
bursts (Pallavicini et al. 1985, Lang & Willson 1986a, 1988a, White & Franciosini 1995, Large 
et al. 1989). 



12 



Giidel 



6 RADIO FLARES AND CORONAL HEATING 

6.1 Is Quiescent Emission Composed of Flares? 

The question of the nature of stellar quiescent radio emission has defined one of the most 
fascinating aspects of stellar radio astronomy. A number of observations seem to suggest that 
flares play an important role - quiescent emission could simply be made up of unresolved flares. 
Inevitably, this question relates to the physics of coronal heating and the presence of X-ray 
coronae. 

Very Long Baseline Interferometry (VLBI) studies of RS CVn binaries suggest that "flare 
cores" progressively expand into a large-scale magnetosphere around the star, radiating for sev- 
eral days ("flare remnant emission"). The electron distributions then evolve from initial power- 
law distributions as they are subject to collisional losses affecting the low-energy electrons and 
to synchrotron losses affecting the high-energy electrons. Time-dependent calculations predict 
rather flat spectra similar to those observed (Chiuderi Drago & Franciosini 1993, Franciosini & 
Chiuderi Drago 1995, Torricelli-Ciamponi et al. 1998). 

6.2 Microflaring at Radio Wavelengths 

Very long quiescent episodes with little flux changes impose challenges for flare-decay models. 
Limitations on the collisional losses (equation 25) require a very low ambient electron density 
(Massi & Chiuderi-Drago 1992). Frequent electron injection at many sites may be an alternative. 
Based on spectral observations, White & Franciosini (1995) suggest that the emission around 
1.4 GHz is composed of a steady, weakly polarized broad-band gyrosynchrotron component plus 
superimposed, strongly and oppositely polarized, fluctuating plasma emission that is perceived 
as quasi-steady but that may occasionally evolve into strong, polarized flare emission (Simon et 
al. 1985, Lestrade et al. 1988, Jones et al. 1996, Dempsey et al. 2001; Figure 4). 

RS CVn and Algol-type binaries also reveal significant, continuous gyrosynchrotron variability 
on time scales of ~10— 90 minutes during ~30% of the time, with an increasing number of events 
with decreasing flux (Lefevre et al. 1994). These flare- like events may constitute a large part of 
the quiescent gyrosynchrotron emission. However, "microbursts" or "nanoflares" with durations 
on the order of seconds to a few minutes are usually not detected at available sensitivities (Kundu 
et al. 1988, Rucinski et al. 1993) although Willson & Lang (1987) found some rapid variability 
on time scales between 30 s and 1 h that they interpreted as being due to variable absorption 
by thermal plasma. 

6.3 Radio Flares and the Neupert Effect 

The chromospheric evaporation scenario devised for many solar flares assumes that accelerated 
coronal electrons precipitate into the chromosphere where they lose their kinetic energy by colli- 
sions, thereby heating the cool plasma to coronal flare temperatures and evaporating it into the 
corona. The radio gyrosynchrotron emission from the accelerated electrons is roughly propor- 
tional to the injection rate of electrons, whereas the X-ray luminosity is roughly proportional to 
the accumulated energy in the hot plasma. To first order, one expects 

L R {t) oc j t Lx(t), (27) 

a relation that is known as the "Neupert Effect" (Neupert 1968) and that has been well observed 
on the Sun in most impulsive and many gradual flares (Dennis & Zarro 1993). The search for 
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Figure 4: Light curves of HR 1099 (left) and UX Ari (right) obtained in the two senses of circular 
polarization. The brighter of the two polarized fluxes varies rapidly and has been interpreted 
as 100% polarized coherent emission superimposed on a gradually changing gyrosynchrotron 
component (White & Franciosini 1995; figures from SM White). 



stellar equivalents has been a story of contradictions if not desperation (Figure 5). A first 
breakthrough came with simultaneous EUV (a proxy for X-rays) and optical observations (a 
proxy for the radio emission) of a flare on AD Leo (Hawley et al. 1995) and radio + X-ray 
observations of flares on UV Cet (Giidel et al. 1996a; Figure 5a). The latter revealed a relative 
timing between the emissions that is very similar to solar flares. Also, the energy ratios seen in 
non-thermal and thermal emissions are similar to solar analogs but, perhaps more interesting, 
they are also similar to the corresponding ratio for the quiescent losses. These observations 
suggest that high-energy particles are the ultimate cause for heating through chromospheric 
evaporation not only in flares, but also in the "quiescent" state. In retrospect, further suggestive 
examples of the Neupert effect can be found in the published literature, most notably in Vilhu 
et al. (1988), Stern et al. (1992), Brown et al. (1998), or Ayres et al. (2001). 

It is important to note that the Neupert effect is observed neither in each solar flare (50% of 
solar gradual flares show a different behavior; Dennis & Zarro 1993), nor in each stellar flare. 
Stellar counter-examples include an impulsive optical flare with following gradual radio emission 
(van den Oord et al. 1996), gyrosynchrotron emission that peaks after the soft X-rays (Osten 
et al. 2000), and an X-ray depression during strong radio flaring (Giidel et al. 1998). Note also 
that complete absence of correlated flaring has been observed at radio and UV wavelengths (e.g., 
Lang & Willson 1988b). 
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Figure 5: Left: Neupert effect seen in an M dwarf star, compared with a solar example in the 
upper panel (Giidel et al. 1996a). Right: Neupert effect seen in an RS CVn binary during a 
large flare (Giidel et al. 2002). 

6.4 The Correlation between Quiescent Radio and X-Ray Emissions 

Radio detections are preferentially made among the X-ray brightest stars (White et al. 1989a), 
a result that is corroborated by new, unbiased all-sky surveys (Helfand et al. 1999). For RS 
CVn and Algol binaries, one finds a correlation between the radio and X-ray luminosities, Lr oc 
l lo-i.3 ( Drake et al 1986; 1989; 19 g 2) Fox et a l. 1994, Umana et al. 1998), whereas for late-type 

dwarfs, a rather tight linear correlation appears to apply (Giidel et al. 1993). In fact, several 
classes of active stars and solar flares follow a similar relation. Overall, for 5—8 GHz emission, 

Lx/L R « 10 15±1 [Hz] (28) 

(Giidel & Benz 1993, Benz &; Giidel 1994; Figure 6). Phenomenologically, the relation simply 
suggests that both radio and X-ray emissions are activity indicators that reflect the level of 
magnetic activity. The most straightforward physical model, proposed by Drake et al. (1989, 
1992), assumes that hot plasma emits both thermal X-rays and non-thermal gyrosynchrotron 
radiation, but this model predicts steep (oc optically thin spectra that are not observed. 

Chiuderi Drago & Franciosini (1993) argued that the build-up of magnetic energy in magnetic 
loops is oc B 2 and proportional to the gradient of the turbulent velocity (assumed to be similar in 
all stars). Since most of the energy in the magnetic fields is eventually radiated away in X-rays, 
Lx oc B 2 V , where V is the source volume. A power-law electron population in the same loop 
emits optically thin radio radiation as Lr oc f?-°- 22 +°- 9<5 y for given v (equation 17). Therefore, 
L R oc l^ 455 - - 11 ^ 1 - 11 - - 45 ' 5 . For 5 = 3, the dependence on V is small and L R oc L^ 24 , close to 
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Figure 6: Correlation between quiescent radio and X-ray luminosities of magnetically active 
stars (symbols) and solar flares (letters; after Benz & Giidel 1994). 

the reported relations. 

If the quiescent emission is made up of numerous solar-like coronal flares, then the high-energy 
particles carry their energy downward to heat chromospheric plasma. For a steady-state situation 
in which the electron injection rate is balanced by an unspecified electron loss mechanism, one 
derives a relation between the synchrotron losses, the particle life time, and the energy losses in 
X-rays. Using the empirical correlation, Giidel & Benz (1993) estimated an electron life time 
that indeed corresponds quite well to the observed radio variability time scales. Holman (1985, 
1986), and Airapetian & Holman (1998) presented detailed model calculations for the alternative 
situation of simultaneous heating and particle acceleration in current sheets. 

7 STELLAR MAGNETIC FIELDS 

Stellar radio astronomy has provided invaluable information on stellar coronal magnetic fields 
not accessible by any other methods. We briefly describe a few principal results. Electron 
cyclotron maser emission can be used to estimate the field strength in the source region, namely 
in converging, lower-coronal fields. Because the radiation is emitted at the fundamental (s = 1) 
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or the second harmonic (s = 2) of the gyrofrequency, B = i/h z /(s x 2.8 x 10 6 ) pa 250 — 500 G 
(Giidel et al. 1989a, Bastian et al. 1990). 

Coronal gyroresonance emission detected at 8—23 GHz is optically thick at harmonics s = 3—5 
of vb (see "Gyromagnetic Emission Mechanisms," above), hence B = ^h z /(s x 2.8 x 10 6 ) pa 
600 — 2070 G in low-lying loops in active regions (Giidel <fe Benz 1989). 

Large coronal sizes as seen by VLBI restrict field strengths to typically 10—200 G if the 
approximate dipolar extrapolation to the photosphere should not produce excessive surface field 
strengths (Benz et al. 1998, Mutel et al. 1998, Beasley & Giidel 2000). 

Very high synchrotron brightness temperatures Tb constrain, together with the power-law 
index S, the characteristic harmonic numbers s and therefore the field strength (Dulk & Marsh 
1982). Lestrade et al. (1984a, 1988), Mutel et al. (1984, 1985), and Slee et al. (1986) found B 
between pa 5 G and several tens of G for RS CVn halo sources, and a few tens to several hundred 
G for the core. Further, the electron energy is e = kTf, for optically thick emission, and e > kTt> 
for optically thin emission. Equation 8 for synchrotron emission (s 2> 1) then implies 

B < 2lT1 f eC ]" ^ Pa 3.6 x 10~ 7 - — ^ [G] (29) 

~ e(m e c 2 + kT b ) 2 (1 + 1.7 x 10- 10 T 6 ) 2 L J v ; 

where the < sign applies for optically thin emission. Alternatively, from Dulk & Marsh (1982) 
one finds for gyrosynchrotron emission T& < T c g- = io 6 ' 15 ~ ' 85<5 (^/f?) ' 5+0 ' 085<5 . Such arguments 
lead to field strengths of a few tens to a few hundred G in RS CVn binaries (Mutel et al. 1985, 
Lestrade et al. 1988). 

The synchrotron turnover frequency f pca k determines B through B = 2.9 x W 13 Vp Cak 9 4 S~ 2 
where 6 is the source diameter in arcsec, S is the radio flux density in mJy at f pea k which is 
given in units of GHz here (Lang 1999, after Slish 1963). For a typical RS CVn or dMe coronal 
source with a size of 1 milliarcsec (mas) as measured by VLBI, a turnover at ^ pea k ~ 5 GHz and 
a flux of 10 mJy, magnetic fields of up to a few hundred G are inferred. For dipolar active region 
gyrosynchrotron models, White et al. (1989b) derived B pa 150fp Cak (f pca k m GHz), which again 
implies field strengths of several hundred G for turnover frequencies in the GHz range. 

Full spectral modeling of the magnetic field strength B, the non-thermal particle density, 
the geometric size, and the electron power-law index S can constrain some of these parameters. 
Umana et al. (1993, 1999) and Mutel et al. (1998) found B pa 10 - 200 G in their models for 
core plus halo structures in Algol- type stars. 



8 RADIO CORONAL STRUCTURE 



8. 1 Introduction 

Although the photospheric filling factor of kG magnetic fields on the Sun is small, it can exceed 
50% on late-type active stars. On the Sun, magnetic fields can rapidly expand above the tran- 
sition region and thus drop below kG levels. On magnetically active stars, such divergence is 
prevented by adjacent magnetic flux lines so that strong magnetic fields may exist also at coro- 
nal levels. The effective scale height of the coronal magnetic field however also largely depends 
on the structure of the magnetic flux lines, i.e., on whether they are compact loops with short 
baselines, long loops connecting distant active regions, or large dipolar magnetospheres anchored 
at the stellar poles. This issue is unresolved. Arguments for small-scale coronal active regions 
as well as for star-sized global magnetospheres have been put forward (see discussions in White 
et al. 1989b, Morris et al. 1990, and Storey & Hewitt 1996). 
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8.2 Radio Eclipses and Rotational Modulation 

Eclipses and rotational modulation offer reliable information on radio source geometries, but 
neither are frequently seen. The radio sources may be much larger than the eclipsing star, or 
they may not be within the eclipse zone. Complete absence of radio eclipses has, for example, 
been reported for AR Lac (Doiron & Mutel 1984), Algol (van den Oord et al. 1989, Mutel et 
al. 1998), and YY Gem (Alef et al. 1997). Positive detections include the Algol system V505 
Sgr (Gunn et al. 1999) and the pre-cataclysmic system V471 Tau (Patterson et al. 1993, Lim 
et al. 1996c). The former surprisingly shows both a primary and a secondary eclipse although 
one of the components is supposed to be inactive. The authors suggested that a radio coronal 
component is located between the two stars, a conjecture that gives rise to interesting theoretical 
models such as diffusion of magnetic fields from the active to the inactive star (on unreasonably 
long time scales, however), or radio emission from the optically thin mass accretion stream 
(which is found to be too weak), or field shearing by the inactive companion. The radio emission 
of V471 Tau may originate from magnetic loops that extend to the white dwarf where they 
interact with its magnetic field (Lim et al. 1996c). 

Radio rotational modulation is often masked by intrinsic variability (Rucinski 1992). Well- 
documented examples are the RS CVn binaries CF Tuc (Budding et al. 1999) for which Gunn et 
al. (1997) suggested the presence of material in the intrabinary region, and the RS CVn binary 
UX Ari, which appears to have radio-emitting material concentrated above magnetic spots on 
the hemisphere of the K subgiant that is invisible from the companion (Neidhofer et al. 1993, 
Elias et al. 1995, Trigilio et al. 1998). The clearest main-sequence example is AB Dor, with two 
emission peaks per rotation seen repeatedly over long time intervals. Both peaks probably relate 
to preferred longitudes of active regions (Slee et al. 1986, Lim et al. 1992, Vilhu et al. 1993). The 
shape of the maxima suggests that the radiation is intrinsically directed in the source, which is 
plausible for synchrotron emission by ultrarelativistic electrons (Lim et al. 1994). 

8.3 Very Long Baseline Inter ferometry 

Interferometry at large baselines has long been a privilege of radio astronomy and has proven 
versatile for numerous applications (Mutel 1996). Various intercontinental Very Long Baseline 
Interferometry (VLBI) networks have been arranged, and dedicated networks such as the US 
Very Large Baseline Array (VLBA) or the UK MERLIN network are routinely available. Apart 
from astrometric applications with sub-mas accuracy (Lestrade et al. 1990, 1993, 1995, 1999, 
Guirado et al. 1997), VLBI provides coronal mapping at ^ 1 mas angular resolution. 

One of the principal, early VLBI results for RS CVn and Algol-like binaries relates to evidence 
for a compact core plus extended halo radio structure of a total size similar to the binary system 
size (Mutel et al. 1984, 1985, Massi et al. 1988, Trigilio et al. 2001). Compact cores appear 
to be flaring sites, whereas halo emission corresponds to quiescent, low-level radiation, perhaps 
from decaying electrons from previous flares. There is evidence for non-concentric, moving, or 
expanding sources (Lestrade et al. 1988, Trigilio et al. 1993, Lebach et al. 1999, Franciosini et 
al. 1999, Lestrade et al. 1999). 

Pre-main sequence stars are attractive VLBI targets. Although some sources are overresolved 
at VLBI baselines (Phillips et al. 1993), others have been recognized as components in close 
binaries such as HD 283447 (Phillips et al. 1996); but most of the nearby weak- lined T Tau 
stars (see "Evolution to the Main Sequence," below) reveal radio structures as large as 10i?*, 
indicative of extended, global, probably dipolar-like magnetospheres somewhat similar to those 
seen in RS CVn binaries (Phillips et al. 1991, Andre et al. 1992). Perhaps also related to youth, 
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Figure 7: (a) VLBA image of the dMe star UV Cet; the two radio lobes are separated by 
about 1.4 mas, while the best angular resolution reaches 0.7 mas. The straight line shows the 
orientation of the putative rotation axis, assumed to be parallel to the axis of the orbit of UV Cet 
around the nearby Gl 65 A. The small circle gives the photospheric diameter to size, although 
the precise position is unknown (after Benz et al 1998). (b) VLBA image of the Algol binary. 
The most likely configuration of the binary components is also drawn. The radio lobes show 
opposite polarity (Mutel et al 1998; reproduced with permission of the AAS.) 

magnetospheres around Bp/Ap stars have also been observed with VLBI and have strongly 
supported global magnetospheric models (Phillips & Lestrade 1988, Andre et al. 1991; see "Stars 
at the Interface Between Hot Winds and Coronae," below). 

VLBI techniques have been more demanding for single late-type dwarf stars, owing both to 
lower flux levels and apparently smaller coronal sizes. Some observations with mas angular 
resolutions show unresolved quiescent or flaring sources, thus constraining the brightness tem- 
perature to Tfe > 10 10 K (Benz & Alef 1991, Benz et al. 1995), whereas others show evidence 
for extended coronae with coronal sizes up to several times the stellar size (Alef et al. 1997, 
Pestalozzi et al. 2000). 

The dMe star UV Cet is surrounded by a pair of giant synchrotron lobes, with sizes up 
to 2.4 x 10 10 cm and a separation of 4—5 stellar radii along the putative rotation axis of the 
star, suggesting very extended magnetic structures above the magnetic poles (Benz et al. 1998; 
Figure 7a). VLBA imaging and polarimetry of Algol reveals a similar picture with two oppositely 
polarized radio lobes separated along a line perpendicular to the orbital plane by more than the 
diameter of the K star (Mutel et al. 1998; Figures 7b, 8b). Global polarization structure is 
also suggested in UX Ari, supporting the view that the magnetic fields are large-scale and well 
ordered (Beasley & Gudel 2000). 
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Figure 8: (a) Equatorial model for the magnetosphere of the young B star SI in p Oph (Andre et 
al. 1988). (b) Sketch for radio emission from a global dipole consistent with the VLBI observation 
shown in Figure 7b (Mutel et al. 1998; reproduced with permission of the AAS.) 

8.4 Magneto spheric Models 

Quite detailed geometric models of large magnetospheres around RS CVn and Algol binaries, T 
Tau stars, and magnetic Bp/Ap stars have been designed based on VLBI results, radio spectra 
and polarization. Common to all is a global, dipole-like structure somewhat resembling the 
Earth's Van Allen belts (Figure 8a). Stellar winds escaping along magnetic fields draw the 
field lines into a current sheet configuration in the equatorial plane. Particles accelerated in 
that region travel back and are trapped in the dipolar-like, equatorial magnetospheric cavity. 
Variants of this radiation belt model, partly based on theoretical work of Havnes & Goertz 
(1984), have been applied to RS CVn binaries (Slee et al. 1987b, Morris et al. 1990, Jones et al. 
1994, Storey & Hewitt 1996), in an optically thick version to Bp/Ap stars (Drake et al. 1987b, 
Linsky et al. 1992) and in an optically thin version to a young B star (Andre et al. 1988). 

Such models are particularly well supported by polarization measurements in RS CVn bina- 
ries. For a given system, the polarization degree p at any frequency is anticorrelated with the 
luminosity, while the sense of polarization changes between lower and higher frequencies. For a 
stellar sample, p is inversely correlated with the stellar inclination angle such that low-inclination 
("pole-on") systems show the strongest polarization degrees (Mutel et al. 1987, 1998, Morris et 
al. 1990). 

When flares are in progress, the two-component core plus halo model appears to correctly 
describe radio spectral properties. Stronger magnetic fields of 80—200 G are inferred for the 
core and weaker fields of 10—30 G for the halo (Umana et al. 1993, 1999, Mutel et al. 1998, 
Trigilio et al. 2001). The frequency-dependent optical depth makes the source small at high 
frequencies (size « i?*, above 10 GHz) and large at small frequencies (size comparable to the 
binary system size at 1.4 GHz; Klein & Chiuderi-Drago 1987, Jones et al. 1995). This effect 
correctly explains the relatively flat, optically thick radio spectra seen during flares. 
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9 STELLAR ROTATION, BINARITY, AND MAGNETIC CYCLES 

9. 1 Radio Emission and Stellar Rotation 

There is little doubt that rotation is largely responsible for the level of magnetic activity in 
cool stars. It is clearly the young, rapidly rotating stars or spun- up evolved stars in tidally 
interacting binaries that define the vast majority of radio-emitting cool stars detected to date. 
The rather peculiar class of FK Com stars, single giants with unusually short rotation periods 
of only a few days, are among the most vigorous emitters of gyrosynchrotron emission (Hughes 
& McLean 1987, Drake et al. 1990, Rucinski 1991, Skinner 1991), including extremely luminous 
and long flares (Slee et al. 1987b, Bunton et al. 1989). However, large surveys of RS CVn binary 
systems find at best a weak correlation between radio luminosity and rotation parameters (Mutel 
& Lestrade 1985, Morris & Mutel 1988, Caillault et al. 1988, Drake et al. 1989, 1992), and it 
may even depend on the optical luminosity class considered. Although Stewart et al. (1988) 
reported a trend of the form Lr oc Rl' 5 v 2 ' 5 (v is the equatorial velocity) for different luminosity 
classes, this result was later challenged (Drake et al. 1989) and may be related to using peak 
fluxes and flares. The absence of a correlation may be related to the coronal saturation regime 
known in X-rays (Vilhu & Rucinski 1983, Mutel & Lestrade 1985, White et al. 1989a): almost 
all radio-detected stars emit at the maximum possible and rotation-independent X-ray level - 
hence perhaps also at the maximum possible radio level. There is obviously a large potential for 
stellar radio astronomy in the less exotic regime below saturation! 

Moving to tighter binary systems in which the components are in (near-)contact, we would 
expect magnetic activity to increase or at least to stay constant. This is however not the case. 
The radio emission of such systems is significantly weaker than that of RS CVn binaries or active 
single stars, a trend that is also seen at other wavelengths (Hughes &: McLean 1984, Rucinski 
Sz Seaquist 1988, Vilhu et al. 1988). Possible physical causes include i) a reduced differential 
rotation, hence a weaker dynamo action (Beasley et al. 1993), ii) a shallower convection zone 
in contact binary systems, and iii) influence by the energy transfer between contacting stars 
(Rucinski 1995). 

9.2 Activity Cycles 

Radio polarization may be a telltale indicator for magnetic activity cycles on stars, analogous 
to the solar 11 year cycle (or 22 years, if the magnetic polarity reversal is considered). Long- 
term measurements in the Ca H&K lines (the HK project; Baliunas 1998) indicate the presence 
of activity cycles with durations of several years in low-activity stars, but irregular long-term 
variations on active stars. In any case, one would expect reversals of the dominant sense of 
polarization to accompany any of these quasi-cycles. The contrary is true, to an embarrassing 
level: After decades of monitoring, many active stars show a constant sign of radio polarization 
throughout, both in quiescence and during flares (Gibson 1983, White et al. 1986, Mutel et al. 
1987, Kundu & Shevgaonkar 1988, Jackson et al. 1989, White & Franciosini 1995, Mutel et al. 
1998) with few exceptions (Bastian & Bookbinder 1987, Willson et al. 1988, Lang 1990). This 
suggests the presence of some stable magnetospheric structures or a predominance of strong 
magnetic fields in one polarity. A concerted and ongoing effort to look for polarity reversals is 
negative at the time of this writing (S. White, private communication). 

Perhaps the most convincing case yet reported in favor of a magnetic cycle is HR 1099 in which 
the average radio flux density correlates with the spot coverage, revealing a possible periodicity 
of 15—20 years (Umana et al. 1995). Massi et al. (1998) reported a surprisingly rapid quasi- 
periodicity in microwave activity (P « 56 ± 4 days) and the accompanying sense of polarization 
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Figure 9: (a) Betelgeuse observed at 43 GHz. The radio photosphere of this star is resolved, with 
an angular resolution of 40 mas. (b) The atmospheric temperature of Betelgeuse as a function of 
radius, observed at different frequencies. (From Lim et al. 1998; reproduced with the permission 
of Nature.) 

(P ps 25 days) in the RS CVn binary UX Ari. Whether such rapid oscillations are related to an 
internal magnetic dynamo cycle remains to be shown. 

10 RADIO EMISSION FROM CHROMOSPHERES AND WINDS 

10.1 Lower Atmospheres of Cool Stars 

Stellar transition regions and chromospheres are expected to be radio sources as well (see "Qui- 
escent Emission from Coronae of Cool Stars," above). Among cool main-sequence stars, the 
slightly evolved mid-F star Procyon is the only such source in the solar vicinity detected to date 
(Drake et al. 1993). In the red giant area, however, outer atmospheres fundamentally change 
their characteristics. Chromospheres become as large as several stellar radii. Such sources are 
now well detected at radio wavelengths (Knapp et al. 1995) and spatially resolved (Skinner et al. 
1997), albeit with some surprises. Because their outer atmospheres are optically thick at radio 
wavelengths, spatially resolved observations provide a direct temperature measurement. Reid 
& Menten (1997) inferred optically thick "radio photospheres" at about 2R* but at subphoto- 
spheric/nonchromospheric temperatures. Indeed, contrary to chromospheric UV measurements, 
the radio-derived temperature in the nearby supergiant Betelgeuse is seen to drop systemati- 
cally from optical-photospheric levels outward (Lim et al. 1998; Figure 9). This cool material 
completely dominates the outer atmosphere. The authors suggested that cool, photospheric ma- 
terial is elevated in giant convection cells; dust formation in this environment could then drive 
Betelgeuse 's outflow. 

10.2 Winds from Cool Stars 

Could dwarf stars lose mass by a substantial stellar wind (Doyle & Mathioudakis 1991, Mullan 
et al. 1992)? Such a wind, if (partially) ionized, would be a radio source. The very fact 
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that coronal radio flares are seen at GHz frequencies implies that the extended winds must be 
optically thin down to coronal heights, and this places stringent limits on the mass loss rate. 
Sensitive millimeter measurements and theoretical arguments have constrained such mass loss 
for M dwarfs to & KT 13 M Q yr -1 or & lCr 12 M yr _1 for a 10 4 K or a 10 6 K wind, respectively 
(Lim & White 1996, Lim et al. 1996a, van den Oord & Doyle 1997), in agreement with upper 
limits derived from UV observations (Wood et al. 2001). Stringent upper limits to the radio 
emission of solar analogs also help confine the mass loss history of the Sun or solar-like stars 
(Drake et al. 1993, Gaidos et al. 2000). An upper limit to the mass loss during the Sun's main- 
sequence life of 6% indicates that the Sun was never sufficiently luminous to explain the "Young 
Sun Paradox", i.e., the suggestion that the young Sun was more luminous in early times given 
the apparently much warmer climate on Mars (Gaidos et al. 2000). 

Winds become progressively more important toward the red giants, in particular beyond the 
corona-cool wind dividing line. Ionized-mass losses between 10~ 10 - 1O" 9 M yr 1 arc indicated, 
increasing toward cooler and more luminous stars (Drake & Linsky 1986). However, the surface 
mass loss flux is similar in all stars, in fact also similar to the solar wind mass loss flux. The 
coolest, late-M and C (super)giants still support massive winds, but the much weaker radio 
emission indicates that the ionization fraction drops by at least an order of magnitude compared 
to earlier M giants, i.e., their chromospheres must be cool, and their optical depth at radio 
wavelengths becomes small (Drake et al. 1987a, 1991a, Luttermoser & Brown 1992). The same 
holds for technetium-deficient S and M-S giants, although white dwarf companions of some of 
them may sufficiently ionize the outer atmosphere of the giant to become visible at radio wave- 
lengths (Drake et al. 1991b). This latter mechanism is probably also relevant in the interacting 
"symbiotic binaries" that usually consist of a red giant and a white dwarf companion and that 
are radio sources with considerable luminosities. The white dwarf is sufficiently UV strong to 
ionize part of the cool-star wind (Seaquist et al. 1984, 1993, Taylor & Seaquist 1984, Seaquist 
& Taylor 1990). No appreciable radio emission is, however, detected from the similar class of G 
and K-type Barium stars (showing overabundances of Ba and other nuclear-processed elements), 
although they show evidence of white-dwarf companions (Drake et al. 1987c). 

10.3 Ionized Winds and Synchrotron Emission from Hot Stars 

Because OB and Wolf-Rayet stars shed strong ionized winds, they emit thermal-wind radio 
emission (equations 24). Some of the hot-star radio sources are very large and resolved by 
the VLA because the optically thick surface (equation 23) is located at hundreds of stellar 
radii. Under the assumption of a steady, spherically symmetric wind, wind mass-loss rates 
of M ph 10~ 6 - 10~ 5 M y^ 1 are inferred for O and B stars (Scuderi et al. 1998), and M » 
2-4 x 10~ 5 M & yr" 1 for WR stars (Bieging et al. 1982, Leitherer et al. 1995, 1997). The inferred 
mass-loss rate is closely correlated with the bolometric luminosity, logM = (1.15 ± 0.2)logL + C 
(Scuderi et al. 1998), in good agreement with Ha measurements. In the case of colliding wind 
binaries, the thermal radio emission can be further enhanced by contributions from the wind- 
shock zone (Stevens 1995). The radio emission level drops appreciably toward intermediate 
spectral classes of B, A, and F, probably owing to a steep decrease of the ionized mass loss 
(Drake & Linsky 1989). 

It came as a surprise when several OB stars (Abbott et al. 1984, Bieging et al. 1989, Drake 
1990, Contreras et al. 1996) and WR stars (Becker k White 1985, Caillault et al. 1985, Abbott 
et al. 1986, Churchwell et al. 1992, Chapman et al. 1999) were found to show non-thermal, 
synchrotron-like radio spectra, some of them associated with short-term variability (Persi et al. 
1990). Up to 50% of the WR stars appear to be non-thermal sources, and this fraction is half 
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Figure 10: VLA map of the WR 147 system observed at 3.6 cm. Also shown are calculated model 
curves that follow the shock formed by the two-wind interaction (from Contreras & Rodriguez 
1999; reproduced with the permission of the AAS). 

as large for OB stars (Leitherer et al. 1997, Chapman et al. 1999, Dougherty & Williams 2000). 
Hot stars are not thought to produce magnetic fields via a dynamo. Moreover, because the wind 
is optically thick to radio emission out to hundreds of stellar radii, the non-thermal component 
must originate at such large distances as well. A coronal model should therefore not apply 
although highly variable radio emission and a poor correlation with X-ray behavior may suggest 
that magnetic fields play a role in the structuring of the winds (Waldron et al. 1998). Viable 
alternatives include synchrotron emission from electrons accelerated in shocks of unstable winds 
of single stars (White 1985, Caillault et al. 1985), in colliding-wind shocks in massive binaries 
(Eichler & Usov 1993), and in the interaction zone between the thermal wind and a previously 
ejected shell (Luminous Blue Variable phase; Leitherer et al. 1997). Magnetic fields inferred for 
the photospheric level are of order 1 Gauss (Bieging et al. 1989, Phillips & Titus 1990), and up 
to 50 G in synchrotron envelope models (Miralles et al. 1994). 

Shocks are attractive features as they can accelerate electrons by the first-order Fermi mech- 
anism, with a predicted electron energy power-law index 5 = 2 (Bell 1978). Skinner et al. 
(1999), however, found a significant deviation from this model for a colliding-wind binary but 
successfully interpreted their radio spectra with an absorbed synchrotron spectrum from a quasi- 
monoenergetic electron population. The origin and stability of such a population are unclear. 
The colliding-wind model found strong support when it became evident that most non-thermal 
sources are indeed binaries (Dougherty & Williams 2000) and that the non-thermal sources are 
located between the two stars, separate from the wind (Moran et al. 1989, Churchwell et al. 1992, 
Dougherty et al. 1996, 2000, Williams et al. 1997, Chapman et al. 1999). Convincing evidence 
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for colliding winds is available for WR 146 (Dougherty et al. 2000), Cygnus OB2 No 5 (Contreras 
et al. 1997), and WR 147 that all show a thermal source coincident with the WR star and a 
separate bow-shock shaped non-thermal source in the wind collision zone close to a lower-mass 
companion (Figure 10, Contreras k, Rodriguez 1999). If the stars are in a strongly eccentric 
orbit, the wind-shock zone between two stellar components can enter the thermal-wind radio 
photosphere and become absorbed, thus producing strong long-term variability of the radio flux 
(Williams et al. 1990, 1992, 1994, White & Becker 1995, Setia Gunawan et al. 2000). Further 
long-term synchrotron variability may be caused by long-term modulation in the magnetic field 
along the orbit, whereas short-term (~daily) variability may be due to clumps in the wind that 
arrive at the shock (Setia Gunawan et al. 2000). 

The rather inhomogeneous class of Be stars shows evidence for very large (~ 100i?*) high- 
density outflowing disks that have been probed at radio wavelengths. Steeply increasing radio 
spectra and some flux variability are characteristic, but there is considerable debate on the source 
geometry (Taylor et al. 1987, 1990, Drake 1990, Dougherty et al. 1991). 

There have been a number of interferometric observations of non-thermal sources among hot 
stars. Source sizes of order 100.R*, exceeding the size of the optically thick surface of the thermal 
wind, have been measured. They imply brightness temperatures up to at least 4 x 10 K and 
thus further support the non-thermal interpretation (Felli et al. 1989, 1991, Felli & Massi 1991, 
Phillips & Titus 1990). 

11 STAR FORMATION AND THE SOLAR-STELLAR CONNECTION 

11.1 Radio Emission from Low-Mass Young Stellar Objects 

Large numbers of visible T Tau stars and embedded infrared sources in star forming regions 
are strong radio sources (Garay et al. 1987, Felli et al. 1993). Schematically, pre-main se- 
quence evolution is thought to proceed through four consecutive stages with progressive clearing 
of circumstellar material (Shu et al. 1987): i) Cold condensations of infalling molecular mat- 
ter, forming a hydrostatic low-luminosity protostellar object, with the bulk mass still accreting 
("Class source"), ii) Formation of a deeply embedded protostar ("Class I source") through 
further accretion via a massive accretion disk, associated with strong polar outflows, iii) The 
"classical T Tau" (CTT, "Class II") phase with an optically visible star accompanied by a thick 
accretion disk, a weak outflow, and possibly a weakly ionized wind; and iv) the "weak-lined 
T Tau" (WTT, "Class III") phase at which disk and circumstellar material have largely been 
dissipated, and the star approaches the main sequence. The evolutionary sequence is somewhat 
controversial (see below) and may in fact describe the evolutionary phase of the circumstellar 
material rather than of the star itself. 

At first inspection, and especially in the outer reaches of the molecular clouds such as the 
p Oph cloud (Andre et al. 1987, 1988, Stine et al. 1988, Magazzu et al. 1999), one encounters 
predominantly WTTS that appear to have evolved past the CTT phase. How early in its 
infancy can a star be to develop strong radio signatures? Whereas deeper surveys of p Oph have 
accessed several deeply embedded infrared sources with high radio luminosities (logLR > 15, 
some of which are class I sources, Andre et al. 1987, Brown 1987, Leous et al. 1991, Feigelson et 
al. 1998) many of the radio-strong WTTS in the p Oph dark cloud are significantly closer to the 
sites of current star formation, and therefore younger, than the typical radio-weaker WTTS and 
CTTS. It seems that radio-strong WTTS evolve directly from embedded protostars. Andre et 
al. (1992) speculated that in some cases strong fossil magnetic fields accelerate both dissipation 
of circumstellar material and formation of large magnetospheric structures on short time scales 
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Figure 11: Observation of the T Tau system at 6 cm with the MERLIN interferometer. Right- 
and left-hand circularly polarized components are shown dashed and solid, respectively. The 
offset of the polarized flux centroids in T Tau(S) is interpreted in terms of polarized outflows 
(Ray et al. 1997; reproduced with the permission of Nature.) 



of ~ 10 6 years. 

Genuine, embedded class I protostars have most often been detected as thermal sources, and 
this emission is predominantly due to collimated thermal winds or jets. These jets are probably 
ionized by neutral winds that collide with the ambient medium at distances of around 10 AU 
and that are aligned with molecular outflows (Bieging & Cohen 1985, Snell et al. 1985, Brown 
1987, Curiel et al. 1987, 1989, 1990, Rodriguez et al. 1989, 1990, 1995, Morgan et al. 1990, Marti 
et al. 1993, Garay et al. 1996, Suters et al. 1996, Torrelles et al. 1997). In the case of more 
massive stars, the radio emission can also originate from optically thick or thin compact HII 
regions (e.g., Hughes 1988, Estalella et al. 1991, Gomez et al. 2000), or from ionized winds (Felli 
et al. 1998); even (thermal or non-thermal) gyrosynchrotron emission has been proposed given 
the high brightness temperatures, small sizes, variability, and negative spectral indices of some 
sources (Hughes 1991, Hughes et al. 1995, Garay et al. 1996). For a review of thermal radio 
jets driving outflows and Herbig-Haro objects, see, e.g., Anglada (1995, 1996), and Anglada et 
al. (1998). Ionized circumstellar material and winds easily become optically thick and therefore 
occult any non-thermal, magnetic emission from close to the star. However, the discovery of 
polarization in T Tau(S) (Phillips et al. 1993), in IRS 5 (Feigelson et al. 1998), in protostellar 
jet sources (Yusef-Zadeh et al. 1990) and the jet outflows themselves (Curiel et al. 1993, Hughes 
1997, Ray et al. 1997; Figure 11), as well as variability and negative spectral indices in T Tau(S) 
(Skinner h Brown 1994) provided definitive evidence for magnetic fields and particle acceleration 
in these class I objects. 

Moving toward the youngest accreting or class sources, centimetric thermal radio detections 
probably again relate to jets/collimated winds that drive massive outflows, whereas the bulk of 
the emitted power leaves the system at mm/submm wavelengths from cold ( <^ 20 K) dust, a 
defining property of class protostars (Andre et al. 1993). Several objects have been detected, 
with ages of the order of only 10 4 years and central stellar masses of only pa 0.05M Q (Leous et 
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al. 1991, Andre et al. 1993, Bontemps et al. 1995, Yun et al. 1996, Andre et al. 1999, Gibb 1999, 
Reipurth et al. 1999), marking the very beginning of the protostellar accretion phase. Radio 
emission is thus a sensitive tracer for the presence of an embedded, nascent but already formed 
protostellar core as opposed to a contracting cloud fragment (see, e.g., Yun et al. 1996). 

11.2 Evolution to the Main Sequence 

Early VLA surveys of CTTS quickly recognized their thermal wind-type emission with rising 
spectra and large angular sizes (Cohen et al. 1982, Bieging et al. 1984, Cohen & Bieging 1986, 
Schwartz et al. 1984, 1986). The enormous kinetic wind energy derived under the assumption 
of a uniform spherical wind suggests anisotropic outflows while structural changes in the radio 
sources indicate variable outflows, probably along jet-like features (Cohen &: Bieging 1986). 
Mass loss rates are estimated to be ^ 1O _7 M yr -1 (Andre et al. 1987). It is important to note 
that the thermal radio emission says nothing about the presence or absence of stellar magnetic 
fields. CTTS do show a number of magnetically induced phenomena, but whatever the possible 
accompanying radio emission, it is thought to be absorbed by the circumstellar ionized wind, 
unless huge magnetospheric structures reach beyond the optically thick wind surface (Stine et 
al. 1988, Andre et al. 1992). 

The CTT stage is pivotal for planetary system formation as massive accretion disks are 
present. Their dust emission dominates the systemic luminosity at millimeter wavelengths. 
Large (>1000 AU) molecular and dust features have been mapped at such wavelengths; how- 
ever, observations of the solar-system sized inner disks (100 AU) have been challenging and 
require high angular resolution. At cm wavelengths, thermal collimated outflows may become 
dominant over dust disks but both features can be mapped simultaneously in some cases, reveal- 
ing two orthogonal structures, one a jet and one an edge-on disk (Rodriguez et al. 1992, 1994, 
Wilner et al. 1996). 

At ages of typically 1 — 20 x 10 6 years, most T Tau stars dissipate their accretion disks 
and circumstellar material and become similar to main-sequence stars albeit at much higher 
magnetic activity levels, probably induced by their short rotation periods (O'Neal et al. 1990, 
White et al. 1992a). The presence of huge flares (Feigelson &: Montmerle 1985, Stine et al. 

1988, Stine & O'Neal 1998), longer-term variability, and falling spectra clearly point to non- 
thermal gyrosynchrotron emission (Bieging et al. 1984, Kutner et al. 1986, Bieging & Cohen 

1989, White et al. 1992a, Felli et al. 1993, Phillips et al. 1996) analogous to more evolved 
active stars. Conclusive radio evidence for the presence of solar-like magnetic fields in WTTS 
came with the detection of weak circular polarization during flare-like modulations, but also in 
quiescence (White et al. 1992b, Andre et al. 1992, Skinner 1993). Zeeman measurements confirm 
the presence of kGauss magnetic fields on the surface of some of these stars (e.g., Johns-Krull 
et al. 1999). Extreme particle energies radiating synchrotron emission may be involved, giving 
rise to linear polarization in flares on the WTT star HD 283447 (Phillips et al. 1996). VLBI 
observations showing large (~10i?*) magnetospheric structures with brightness temperatures of 
T b ss 10 9 K fully support the non-thermal picture (Phillips et al. 1991). 

As a WTT star ages, its radio emission drops rapidly on time scales of a few Million years 
from luminosities as high as 10 18 erg s _1 Hz _1 to values around or below 10 15 erg s _1 Hz _1 at ages 
beyond 10 Myr. Young age of a star is thus favorable for strong radio emission (O'Neal et al. 

1990, White et al. 1992a, Chiang et al. 1996), whereas toward the subsequent Zero- Age Main- 
Sequence (ZAMS) stage it is only the very rapid rotators that keep producing radio emission at 
a 10 15 erg s^Hz" 1 level (Carkner et al. 1997, Magazzu et al. 1999, Mamajek et al. 1999). 
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11.3 Low- Mass Stars on the Main Sequence 

The pace at which radio emission decays toward the ZAMS has prevented systematic radio 
detections of cool main-sequence stars beyond 10—20 pc. Much of the evolutionary trends 
known to date have been derived from fairly small and very select samples of extraordinary stars. 
Although many nearby M dwarfs are probably quite young and located near the ZAMS, their 
ages are often not well known. An interesting exception is the proper-motion dM4e companion 
to AB Dor, Rstl37B, detected as a surprisingly luminous steady and flaring radio star (Lim 
1993, Beasley & Cram 1993). With an age of ~ 5 - 8 x 10 7 years, it may indicate that stars 
approaching the main sequence go through a regime of strongly enhanced magnetic activity. The 
AB Dor pair is a member of the Local Association, also known as the Pleiades Moving Group, 
a star stream of an age (~50— 100 Million years) corresponding to near-ZAMS for F— K stars. 
A handful of further stream members with vigorous radio emission are now known in the solar 
neighborhood, most notably PZ Tel (cited in Lim & White 1995), HD 197890 (Kl V; Robinson 
et al. 1994), HD 82558 (Kl V, Drake et al. 1990), EK Dra (GO V, Giidel et al. 1994, 1995a), 
and 47 Cas (F0 V + G V; Giidel et al. 1995c, 1998). 

These objects are analogs of young open cluster stars with the observational advantage of 
being much closer. Clusters, however, are preferred when more precise stellar ages or large 
statistics are required. Some young clusters house a select group of ultra-fast rotators, stars at 
the extreme of dynamo operation with rotation periods ^ 1 day. The absence of any detections 
in Bastian et al.'s (1988) Pleiades radio survey suggests that their flare properties are similar to 
solar neighborhood stars rather than to the much more energetic outbursts occasionally seen in 
star forming regions. Flaring or quiescent radio emissions have been detected from G— K-type 
members of both the Pleiades (Lim & White 1995) and the Hyades (White et al. 1993, Giidel 
et al. 1996b) although these examples, in part ultra-fast rotators at saturated activity levels 
and in part binary systems, clearly represent only the extreme upper end of magnetic activity. 
Radio emission of normal, single solar analogs rapidly declines to undetectable levels after a few 
hundred Million years (Gaidos et al. 2000). 

11. 4 Herbig Ae/Be stars 

The evolution of intermediate- mass (3 — 2OM ) stars is quite different from, and much faster 
than, that of low-mass stars as they still accrete while already on the main sequence. Given 
their intermediate spectral range, it is of great interest to know whether Herbig Ae/Be pre- main 
sequence stars support convective, magnetic dynamos or whether they resemble more massive 
wind sources. Radio emission is the ideal discriminator. A wind-mass loss interpretation is 
compatible with the expected mass loss rates of 10~ 6 - lO^M^r" 1 (Giidel et al. 1989b, Skinner 
et al. 1993). This interpretation is supported by the large radio sizes (order of 1") and the 
absence of circular polarization or strong variability. The radio luminosity is also correlated 
with the stellar temperature and bolometric luminosity. This is expected because wind mass- 
loss rates increase toward higher stellar masses (Skinner et al. 1993). New radio observations 
complemented by millimeter measurements further indicate the presence of substantial dust 
envelopes (Di Francesco et al. 1997). 

Non-thermal gyrosynchrotron sources exist as well among Herbig stars, although they are the 
exception (Skinner et al. 1993). The evidence comes primarily from negative spectral indices, 
including the extremely X-ray strong proto-Herbig star EC95 = S68-2 that further supports a 
coronal model based on its Lx/Lr ratio (Smith et al. 1999). 
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12 STARS AT THE INTERFACE BETWEEN HOT WINDS AND CORO- 
NAE 

Owing to missing convection, no dynamo-generated magnetic fields are expected on stars earlier 
than spectral type F, nor should there be massive ionized winds on main-sequence stars of 
spectral type B and later, given the weak radiation pressure. Indeed, few detections have been 
reported as early as spectral type FO, although some have been found at quite high luminosities 
compatible with the gyrosynchrotron mechanism, including main-sequence candidates (Giidel et 
al. 1995b) and the supergiant a Car (Slee & Budding 1995). Brown et al. (1990) surveyed a 
number of normal A-type stars and found stringent upper limits to any radio emission and thus 
to the mass-loss rate - very much in agreement with expectations. However, as it turns out, this 
spectral range is shared by some of the more provocative radio detections. 

Many "magnetic chemically peculiar" Bp/Ap stars maintain considerable radio emission (logL^ 
« 15 — 18; Drake et al. 1987b), including very young objects such as SI in p Oph (Andre et al. 
1988). They all relate to the hotter (09— AO spectral type) He-strong and He-weak/Si-strong 
classes, whereas the cooler (A type) SrCrEu peculiarity classes remain undetected (Drake et al. 
1987b, Willson et al. 1988, Linsky et al. 1992, Leone et al. 1994). Similarly, "non-magnetic" Am 
and HgMn stars remain undetected despite recent claims that these stars may have magnetic 
fields as well (Drake et al. 1994). The emission mechanism is likely to be gyrosynchrotron as 
judged from flat spectra (Linsky et al. 1992, Leone et al. 1996), high brightness temperature 
(Phillips & Lestrade 1988), variability, and sometimes moderate polarization (Linsky et al. 1992). 

The presence of kilogauss magnetic fields on magnetic chemically peculiar stars has been known 
since 1947 (Babcock 1947). The fields are generally assumed to be of global, dipolar topology. 
There is little photospheric motion that could stir magnetic footpoints, but weak winds could 
draw the magnetic field lines into an equatorial current sheet, thus producing a global "van Allen 
Belt" magnetospheric structure as described in "Magnetospheric Models," above. Estimated 
non-thermal source radii are a few stellar radii, confirmed by VLBI observations (Phillips & 
Lestrade 1988, Andre et al. 1991) that also conclusively established the non-thermal nature of 
the radio emission, with brightness temperatures of TJ, 10 8 — 10 9 K. 

The wind-controlled magnetospheric model is further supported qualitatively by a parameter 
dependence of the form Lr oc M 0,38 B^s ^rot 32 where M is the estimated wind mass-loss rate 
for the spectral type, B T111S is the root-mean square value of the longitudinal magnetic field, 
and P ro t is the stellar rotation period (the dependence on the latter is marginal; Linsky et al. 
1992). Evidence for rotational modulation probably due to field misalignment has been found 
by Leone (1991), Leone & Umana (1993), and Lim et al. (1996b), the latter authors reporting 
modulation of the polarization degree and sign. A surprise detection was the phase-dependent 
100% polarized radio emission from a Bp star that suggests a strongly beamed, continuously 
radiating electron cyclotron maser (Trigilio et al. 2000). 

A binary class in this spectral range that is thought to be intermediate in evolution between the 
young semi-detached B-type (3 Lyr system that shows radio evidence for large systemic wind mass 
loss (Umana et al. 2000), and evolved normal Algol binaries with non-thermal coronal emission 
is defined by the rather inhomogeneous sample of Serpentid stars or "Peculiar Emission Line 
Algols" (PELAs). These typically consist of an A— B type primary and a F— K type companion 
with strong mass transfer into a geometrically thick accretion disk around the early-type star, 
and a common, thin envelope. They have shown an appreciable level of radio emission, first 
thought to be gyrosynchrotron emission based on their large luminosities (Elias 1990) but later 
suspected to be wind sources given their spectral indices close to the standard wind law (Elias 
& Giidel 1993). 
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13 RADIO EMISSION FROM BROWN DWARFS 

Brown dwarf stars have masses below the critical mass of ss 0.08M Q required for stable hydrogen 
burning in the stellar core. Detecting radio emission that shows variability, polarization, or 
gyrosynchrotron-like spectra from brown dwarfs would provide very strong evidence for the 
existence of surface magnetic fields. However, after brown dwarfs cease to burn deuterium, the 
internal convection decreases, probably diminishing the generation of magnetic fields (Neuhiiuser 
et al. 1999). A deep 3.6 cm survey with the VLA, concentrating on older targets, did indeed 
fail to produce radio detections (Krishnamurthi et al. 1999). One young candidate brown dwarf 
has been reported with strongly variable radio emission (P. Andre, cited in Wilking et al. 1999 
and Neuhiiuser et al. 1999). But again, nature offered more riddles: The first radio-detected 
bona-fide brown dwarf, LP944-20, is in fact quite old (~500 Million years), but showed both 
flaring and quiescent episodes (Berger et al. 2001). Its radio emission is orders of magnitude 
higher than would be expected from the active-stellar relation (28). In any case, the observed 
flaring radio emission is a telltale signature of magnetic fields, and proves that solar-like coronal 
activity is present in substellar objects - even old ones. 

A promising method for detecting sub-stellar objects with sufficiently strong magnetic fields 
involves the operation of the cyclotron maser. In fact, Jupiter's decametric radiation can be as 
strong as solar coherent bursts, although at somewhat longer wavelengths. Searching nearby 
stars for low-frequency coherent bursts may, by positional analysis, reveal companion planets or 
brown dwarfs. Such a study was performed by Winglee et al. (1986), albeit with null results. 
We note in passing that one of their targets, Gl 411/Lalande 21185, has in the meantime been 
shown by astrometric means to possess at least one planetary companion (Gatewood 1996). 

14 STELLAR ASTROMETRY 

VLBI techniques have been used for astrometric purposes, i.e., measurements of positions, proper 
motions, binary orbital motions, and parallaxes, to an accuracy of fractions of a milliarcsecond 
(Lestrade et al. 1990, 1993, 1995, 1999). While one fundamental purpose of such programs is to 
establish the astrometric link between the radio and the optical coordinate reference frames that 
is eventually of great importance for numerous dynamical studies, invaluable astrophysical spin- 
offs were obtained. Astrometry of the Algol triple system (a 2.9-day period A7m+K0 IV close 
binary with mass exchange, in a wide 680 d orbit around a B8 V star) showed unequivocally 
that the radio emission is related to the magnetically active K subgiant. The most startling 
result (confirming earlier, indirect optical polarization measurements) was that the planes of the 
short and the long orbits are orthogonal to each other (Lestrade et al. 1993, 1999), a challenge 
for stability theories. Further important highlights include i) precise measurements of distances 
through parallaxes, including distances to star forming regions, ii) identification of binarity or 
multiple systems, and iii) rapidly moving ejecta after a stellar flare (see Lestrade et al. 1999 
for details). A dedicated pre-Hipparcos VLBI study of the AB Doradus system solved for the 
accurate distance (moving the star to 2/3 its previously assumed distance and placing it on 
the ZAMS) and for a short-period disturbance attributed to the presence of a low-mass object 
(possibly a brown dwarf) in orbit around this star (Guirado et al. 1997). 

15 SUMMARY AND CONCLUSIONS 

Stellar radio astronomy has matured over the past few decades to a science that is indispensable 
for our understanding of stellar atmospheres. Historical milestones include, among many others, 
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the discovery of steady and flaring non-thermal and polarized emission in cool stars, testify- 
ing to the importance of highly energetic processes; the recognition that these phenomena are 
ubiquitous in many classes of convective-envelope stars; observations of very large, apparently 
stable magnetospheric structures, unlike anything known from the Sun, around various types 
of magnetically active stars such as T Tau stars, Bp stars, dMe stars, or RS CVn binaries; the 
discovery of non-thermal emission produced in (wind-collision) shocks of hot-star atmospheres; 
gyromagnetic and flaring emission from deeply embedded protostellar objects, testifying to the 
importance of magnetic fields back to the earliest moments of a stellar life; and flaring radio 
emission from sub-stellar objects not previously thought to support stellar-like convective outer 
envelopes. Radio methodology has become a standard to estimate magnetic fields in cool stars, 
to determine mass loss in stars with ionized winds, to spatially resolve and map structures at the 
milliarcsecond level, and to simply prove the presence of magnetic fields through polarization 
measurements. 

Far from being an auxiliary science to research at other wavelengths, stellar radio astronomy 
should prepare to address outstanding problems to which it has unique access, although more 
sensitive instruments are needed. Questions of particular interest include: Are there relevant 
high-energy processes and magnetic fields in class protostars? Are accretion processes impor- 
tant for the high-energy mechanisms and the generation of large-scale magnetic fields? Are there 
magnetic fields in hot stars, and what role do they play in the winds? Are brown dwarfs usually 
quiescent radio emitters? Do they maintain stable magnetic fields? What is the structure of 
their coronae? Are there intra-binary magnetic fields in close binary stars? How do large mag- 
netospheres couple to the more compact X-ray coronae? Are quiescent coronae fed by numerous 
(micro-)flares? 
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